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Abstract—  Space  Surveillance  Telescope  (SST)  is  a  Defense  Advanced  Research  Projects  Agency  (DARPA) 
program  designed  to  facilitate  the  detection  of  space  debris  in  earth’s  orbit.  In  order  to  achieve  optimal 
performance,  focusing  of  the  telescope  can  be  conducted  by  retrieving  phase  information  in  the  image  to  determine 
the  amount  of  defocus  and  then  moving  the  mirrors  axial  to  shift  the  focal  point.  One  of  its  unique  features  is  that 
operates  with  a  mechanical  shutter  that’s  speed  restricts  the  telescope  to  collecting  long  exposure  imagery.  Long 
exposure  imagery  (» 10ms )  consequently  averages  the  atmosphere,  which  creates  a  point  spread  function  (PSF) 
which  can  mimic  one  that  contains  fixed  aberrations  such  as  focus  and  spherical  error.  The  average  atmosphere 
masks  the  static  aberrations  of  the  telescope  in  the  image  affecting  the  ability  to  achieve  an  optimal  focus.  This 
paper  will  explore  the  joint  estimation  of  the  focus  and  the  atmospheric  seeing  parameter.  The  Cramer-Rao  lower 
bounds  for  variance  are  computed  to  facilitate  an  understanding  of  the  joint  estimation  problem.  These  bounds  will 
serve  to  demonstrate  how  the  average  atmospheric  transfer  function  makes  sensing  a  focus  error  more  difficult  in  the 
presence  of  noise. 

The  views  expressed  are  those  of  the  author  and  do  not  reflect  the  official  policy  or  position  of  the  Department  of 
Defense  or  the  U.S.  Government. 


1.  INTRODUCTION 

The  Department  of  Defense  recently  fielded  an  f/1  Mersenne-Schmidt  telescope  called  the  Space  Surveillance 
Telescope  (SST)  that  advances  the  United  States  space  situational  awareness  (SSA)  capability.  SSA  directly 
supports  the  stated  US  National  Space  Policy  to  “(p)reserve  the  Space  Environment. .  .the  United  States  shall 
develop,  maintain,  and  use  space  situational  awareness  information  from  commercial,  civil,  and  national  security 
sources  to  detect,  identify,  and  attribute  actions  in  space  that  are  contrary  to  responsible  use  and  the  long-term 
sustainability  of  the  space  environment.”  [1] 

The  critical  technology  that  enables  SST’s  6  deg2  field  of  view  (FOV)  camera  is  its  unique  curved  charged 
coupled  device  (CCD)  shown  in 

Fig.  1.  A  high  speed  shutter  was  developed  for  the  camera  with  a  minimum  exposure  time  of  25ms, 
considerably  longer  than  the  \-\0ms  typically  associated  with  a  short  exposure  image  [2]. 


Sensor  Subsystem  (MIT/LL) 


CCD  and  Wide  Field  of  View  Camera  High  Speed  Shutter 


Fig.  1.  SST’s  6  deg  field  of  view  camera  and  high  speed  mechanical  shutter. 
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Currently,  one  of  the  main  challenges  to  optimizing  SSTs  performance  is  to  reduce  the  point  spread 
function  (PSF)  through  focus  and  alignment.  The  pixels  in  the  CCD  are  15 jum  however,  the  pixels  canbe  being  2x2 
binned  due  to  the  expected  atmospheric  blurring.  Fig.  2  shows  the  variations  of  full  width  half  maximum  (FWHM) 
blur  spot  in  terms  of  30 jum  binned  pixels  over  3  months. 


Fig.  2.  Variations  in  the  Full  Width  Half  Maximum  (FWHM)  of  SST’s  Point  Spread  Function  (PSF)  measured  by 
2x2  binned  15 jum  pixels. 

Ideally,  the  PSF  FWHM  would  be  within  one  30 jum  binned  pixel.  Reducing  the  PSF  is  possible  by  accurately 
determining  the  amount  of  focus  error  (and  other  aberrations)  in  the  image  of  a  calibration  star  then  adjusting  the 
focus  and  alignment  to  reduce  the  blur  spot  size.  However  due  to  the  long  exposures,  to  find  unbiased  estimates  of 
the  telescope  aberrations  from  the  star  image  the  atmospheric  effects  must  be  included  in  the  telescope  model. 

2.  TELESCOPE  MODEL 

The  telescope  model  developed  in  this  section  follows  a  similar  model  employed  for  analysis  of  the 
Hubble  Space  Telescope  [3].  The  telescope  is  considered  a  linear  shift  invariant  system  where  the  impulse  response 
of  the  system  will  be  the  PSF.  The  light  propagating  from  the  distance  point  source  (i.e.  a  star)  is  assumed  to  be 
temporally  incoherent.  Therefore,  according  to  Goodman,  the  image  irradiance  in  the  x  detector  plane,  i  (x),  is  the 
convolution  of  point  source  irradiance,  S(x ),  with  telescopes  PSF  [4].  The  parameters  used  in  the  telescope  model 
are  listed  in  Tab.  1. 


Tab.  1.  Telescope  Parameters 


Model  Parameter 

Value 

Center  Wavelength 

500nm 

Telescope  Pupil/Obscuration  Diameter 

3.5m/1.8m 

Telescope  Effective  Focal  Length 

3.5m 

CCD  Pixel  Pitch 

15pm 

Star  Irradiance  per  Frame 

~10A4  photons 

Background  Irradiance  per  Frame 

300  photons 

Grid  Size 

2A11 

The  pupil  function,  A(w,v),  of  the  telescope  is  defined  by  its  annular  aperture  shown  in  Fig.  3a,  where 
u  and  v  are  coordinates  in  the  pupil  plane.  Wave  front  error  caused  by  defocus  is  introduced  into  the  pupil  function 
using  the  Zernike  polynomial  for  defocus,  [5] 

cpA{u,  v)  =  3.464(w2  +  v2)  - 1.732.  (1) 


Distribution  Statement  “A”  (Approved  for  Public  Release,  Distribution  Unlimited) 


The  amount  of  focus  error  is  captured  by  scaling  cp4(u,  v)  with  a  Zernike  coefficient  for  defocus,  Z4 


error(u ,  v)  =  Z4cp4  ( u ,  v). 


(2) 


An  image  of  cp4(u,  v)  scaled  by  a  25  wave  coefficient  is  shown  in  Fig.  3b.  The  aberrations  are  then  represented  in  the 
pupil  plane,  ,  by  the  generalized  pupil  function 


?>(w1)  =  A(w1)  exp  [  j  •  error  (ux )] . 


The  PSF  is  then  computed  as 


KJm)  = 


j  2  71171^ 


=  J\T(Ul  ) eilnmu'  £ T (u2 ) e~ilnm,h- 

u\ 

=YJA(ui)e'r'"’e 

ui  1 

=  \H(m)\2  =H(m)H(m)\ 


jZ4</>4(ux)  ^j27rmu1  e~jZ4<l>4(u2)g-j2xmu2 


(3) 


(4) 


where,  H  is  the  field  in  the  detector  plane  [4].  The  telescope’s  PSF  with  25  waves  of  defocus  is  shown  in  Fig.  3c 
where  m  is  a  pixel  coordinate  in  the  x  detector  plane.  The  large  amount  of  defocus  causes  the  PSF  to  have  an 
annular  shape. 


Fig.  3.  (a)  Pupil  function  used  to  model  SST  (b)  Zernike  polynomial  for  defocus  with  Z4  =  25  waves  (c)  Telescopes 
PSF  with  Z4  =  25  waves  of  focus  error 


For  a  more  complete  model  telescope  model,  the  effect  of  the  finite  square  pixels,  a  =  15 jum,  is  included  in 
the  PSF  where  the  transfer  function  for  the  pixels  and  telescope  respectively  are  represents  as  the  following  digital 
Fourier  transforms, 


^ pm  (u i)  =  T {rect(ax)},  and  (5) 

3-C„p,  (", )  =  F{hop, (m)}’  therefore  (6) 

Ktescope ("*)  =  T  '  (»l  )  ^ pixel  (M1 )}'  (7) 

Images  observed  by  SST  have  been  measured  to  be  shot  noise  dominated,  so  the  image  data,  d(m ),  is  considered  to 
be  Poisson  and  has  a  mean  value  that  is  equal  to  the  irradiance  of  light  in  that  pixel  [2] 


Distribution  Statement  “A”  (Approved  for  Public  Release,  Distribution  Unlimited) 


£[d(m)]  =  i(m). 


(8) 


The  model  for  the  image  irradiance  centered  on  the  optical  axis  is 


i telescope  ( m )  =  Z  ^S(x)htelescope  (m-x)  +  B. 

X 


(9) 


It  includes  additional  terms  to  account  for  the  background  light,  B ,  and  the  total  photons  emitted  from  the  star  per 
integration  time,  6V  The  joint  distribution  of  the  image  data  is  represented  by  the  Poisson  probability  mass  function, 


rwM=n 


~helescope(m)  ' 

1  telescope 


{rn)d(m] 


d(m)\ 


the  associated  log  likelihood  equation  is 


(10) 


L(Z,)  =  \n\J\P[d{m)} 

=  Z  {-telescope  («)  +  d(m)ln  itelescope  (in)  -  In  d{m) !}. 


(11) 


3.  CRAMER-RAO  LOWER  BOUND  (CRLB)  FOR  VARIANCE 

The  Cramer-Rao  lower  bounds  provides  a  theoretical  lower  limit  of  variance  for  estimates  of  the  Zernike  coefficient 
for  defocus,  Z4  [3].  The  bounds  in  Fig.  5  illustrate  that  standard  deviations  of  Z4  on  the  order  of  1CT1  waves  are 
possible  even  in  the  presence  of  long  exposure  atmosphere.  To  determine  the  CRLB  for  estimates  of  the  Zernike 
coefficient  for  defocus  the  Fisher  information,  /(Z4),  is  computed  via  the  following  calculation  [6] 


J(Za)  =  -E 


d2L{z<) 


dZA 


where,  the  CRLB  for  the  variance  of  Z4  is  defined  as 

var(Z4)  >  J  (Z4)  ' . 

The  first  and  second  derivative  of  the  log  likelihood  function  Eq.  (11)  respectively  are 


gL(Z4)_^f  d(m)9l  ^  dhtelescope  (m) 


dZA 


i tele. scope  (^) 


*1 


dZ. 


and 


1  telescope 

q  ^  ^  h telescope  {m)__d{m)a}_(  dh tetocope  (m) 


dZ. 


V  i telescope  (^) 


5Zf 


telescope 


( m ) 


dZ, 


The  resulting  Fisher  information  is 

j{z4)  =  -e 


d2L{z4) 


dZA 


=z 


3 


k  1 telescope  J 


dKle^m) 
v  dZ4  j 


The  derivative  of  the  PSF  with  respect  defocus  is 


(12) 


(13) 


(14) 

(15) 


(16) 
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In  Eq.  (21)  X  is  the  mean  wavelength,  /  is  the  telescope  focal  length,  and  r0  is  the  atmospheric  seeing  parameter. 
The  total  PSF  is  then  computed  as 

Km  (ln)  =  {Mop,  («1 )  'M («!  ) Mam  ( w, )}.  (22) 

Samples  of  the  three  transfer  functions  are  shown  in  Fig.  4  (a-c)  to  illustrate  how  the  pixels  and  atmosphere  reduce 
special  frequency  content  of  the  diffraction  limited  telescope’s  optical  transfer  function  and  due  to  it  Fourier 
transform  relationship  broaden  the  PSF.  The  x-axis  is  shown  in  terms  of  the  spatial  frequency,  ux ,  divided  by  the 

cutoff  frequency,  w0,  for  the  annular  telescope  pupil  function.  As  the  focus  error  increases  J~Copt  (ux )  begins  to 
increasingly  limit  the  spatial  resolution  of  the  telescope  further  broadening  the  PSF. 
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Fig.  4.  (a)  Telescope  Models  Optical  Transfer  Function  (OTF)  with  Z4  =  0  (b)  15 jum  Pixels  Transfer  Function  (c) 
Atmospheric  Transfer  Function  with  r0  =  8 cm 

By  including  the  effects  of  the  atmosphere  in  the  PSF,  the  image  intensity  model  in  Eq.  (9)  becomes 


L  =  Ya6ADKlal(m  -x)  +  B. 

X 

The  elements  of  the  Fisher  information  matrix, 


(23) 


/  = 


/(Z4)  J{z4,r0 ) 

J(Z4,r0)  J(r0 ) 


(24) 


are  calculated  in  order  to  determine  the  CRLB  for  variance  of  Z4  in  the  presence  of  an  average  atmosphere.  Using 
the  log  likelihood  function  in  Eq.  (11)  and  taking  the  second  derivative  of  Eq.  (14)  w.r.t  r0  &  Z4 


d  L(Z4,r0)  

(d(m)6 1  ^ 

Ahtotai  (m)  d(m)6 1 

f  (m)) 

dZ4dr0 

[  LA™)  'J 

3ZA  C«(2(m> 

l  5Z4  j 

5ro  J 

a2L(Z4,r0)  y/  d  (w)  ■ 6>  >  d%^  (m)  d(m)0l 


dr* 


V  L total 

a2L(Z4,r0)_y,r  d  (m)  6>  Q  \  d\^  (m)  d\m)dA  dh Mal  (in) 
dZ}  1  dZ}  i2(m)\  dZ2 


LA™) 


dr'* 


LJim) 


dn 


,  and 


ro  J 


(25) 

(26) 
(27) 


Because  E^d  (m)]  =  itotal  (m)  the  elements  of  the  Fisher  information  matrix  are 


j(Z4,rn)  =  -E 


d2L(Z4,r0 ) 


J(ro)  =  -E 


dZ4dr0 


Ll(z4) 

dr, 2 


=z 


3 


=z 


V  total  J 


^KA™) 

dz4  , 


Sh,oml(™) 


drn 


0  J 


3 


total  J 

f  n  \f 


dh,oA™) 


j(z4)  =  -e 


d2L(Z4) 


dZ ‘ 


dro  2 


,  and 


=Z 


4 


V  *  total  J 


dK,A™) 


dz. 


■4  J 


(28) 

(29) 

(30) 


Then  the  derivatives  of  the  PSF  are 
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(31) 


^hotal  (W)  _  rr- 1  J  rr 

dZ,  1  J 


dH  H+™H* 


dZA 


dZA 


n_T  n_T 

atm^*'  pixel 


=  J-'{j{-2Im  [f[<fi4(u2)A(u2)eJz^] 
dhmal(m)  _  5.73 


H 


rtamtfpixei  k  and 


drn 


The  CRLB  for  variance  is  computed  by  inverting  the  Fisher  information  matrix 


var  (Z4)  co  v(Z4,r0) 

CO  v(Z4,r0)  var(r0) 


(32) 


The  resulting  CRLB  for  the  standard  deviation,  CRLB  <y  ,  of  Z4  is  plotted  in  Fig.  5  for  cases  with  and  without  an 
average  atmosphere  present.  Atmospheric  turbulence  increases  CRLB  <j  and  as  r0  decrease  the  effect  of  the  atmosphere  on  the 
bound  increase.  In  addition  as  Z4  decrease  the  lower  bound  increases.  Therefore,  estimation  of  Z4  should  become  more 
inaccurate  as  the  amount  of  defocus  decreases.  In  addition,  as  r0  decreases  the  CRLB  <y  uniformly  increases  indicating  that 
more  turbulent  atmospheres  will  increase  the  variance  on  estimates  of  defocus. 


Fig.  5.  Cramer-Rao  Lower  Bound  of  the  standard  deviation,  CRLB  cr,  of  estimates  for  the  defocus  parameter,  Z4. 
CRLB  cr  are  shown  for  cases  with  no  atmosphere  in  the  model  and  increasing  atmospheric  seeing  by  changing  r0. 

4.  STAR  SIMULATIONS 

Stars  were  simulated  as  system  impulses,  S(x),  and  then  the  effects  of  the  atmosphere,  telescope,  defocus, 
pixilation,  background  light  and  star  intensity  were  introduced  using  Eq.  (23).  The  analog  images  of  a  star  with  and 
without  atmospheric  effects  are  shown  in  Fig.  6.  (a)  &  (b).  The  pixilated  images  of  those  same  stars  are  picture  in 
Fig.  6.  (c)  &  (d).  Shot  noise  is  simulated  in  the  stars  using  a  Poisson  random  number  generator  the  training  data, 
d(m ),  evaluate  the  performance  of  the  estimators  describes  in  Sec.  8. 
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Fig.  6.  (a)  Simulated  analogue  image  of  a  star  with  Z4  =  18  waves  and  no  atmosphere,  (b)  Simulated  analog  image 
of  a  star  with  Z4  =  18  waves  and  an  average  atmosphere  where  r0  =  8 cm.  (c)  Simulated  digital  image  of  a  star  with 
Z4  =  18  waves  and  no  atmosphere,  (d)  Simulated  digital  image  of  a  star  with  Z4  =  18  waves  and  an  average 
atmosphere  where  r0  =  8 cm. 


5.  PARAMETER  ESTIMATION 


The  method  of  least  squares  (LS)  estimation  was  used  to  estimate  the  Zernike  coefficient  for  defocus,  Z4  from  the 
simulated  star  data.  The  LS  method  is  used  because  of  computer  precision  challenges  encountered  in  the  maximum 
likelihood  estimation  approach  due  to  the  large  background  noise  levels  inside  the  log-likelihood  function.  In 
addition,  the  LS  method  does  not  require  any  parametric  assumptions  making  the  estimator  more  robust  to  changes 
in  the  noise  statistics  [5].  The  intensity  models  from  Eqs.  (9)  &  (23)  are  used  to  define  the  sum  of  squares 


Qtelescope  ^  1  ( telescope  (m)  -  Z  0Ax)hlelescope  ( x-m)-B 


x  V 


,  and 


Qtotal  =  Z  heal (m)  -  Z  dAX)h,o,al  ( X~m)-B 


(33) 

(34) 


In  order  to  estimate  defocus,  the  unknown  photons  per  image,  0X ,  must  also  be  estimated  from  the  data  by  take  the 
derivative  of  Eq.  (33)  and  (34)  and  setting  them  equal  to  zero  to  get  the  generalized  function 


0i  = 


Zft2(m) 


(35) 


Determining  Z4  without  accounting  for  the  atmosphere  is  found  by  the  single  parameter  estimate 

Z4=argmin(e  e), 

Z4 


(36) 
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or  by  accounting  for  the  atmosphere  with  the  joint  estimator 


=  arg  min  (0, (37) 

Z4’r0 

By  conducting  a  numerical  grid  search  of  realistic  values  for  Z4  the  parameter  estimates  that  the  gives  the  LS 

solutions  are  numerically  determined.  Finding  Z4  for  multiple  image  frames  of  simulated  star  data  the  results  are 
used  to  determine  the  sample  mean  and  variance  for  the  LS  estimator. 

To  produce  the  plot  in  Fig.  7a,  training  data,  d(m ),  is  generated  without  shot  noise  and  with  focus  errors 
ranging  from  3-24  waves  in  order  to  determine  the  estimators  biases.  Estimates  of  defocus  using  Eq.  (36)  are  made 
on  simulated  stars  with  and  without  an  average  atmosphere  present.  The  graph  shows  that  when  the  simulated  star 
data  has  an  average  atmosphere,  the  single  parameter  estimator  has  a  defocus  dependent  bias.  In  contrast,  the  results 
of  joint  estimator,  Eq.  (37),  on  the  same  simulated  star  data  with  an  average  atmosphere  present  does  not  have  a 
significant  bias. 

The  joint  estimator  is  used  to  estimate  the  defocus  from  training  data  containing  shot  noise  with  the  mean 
and  standard  deviation  plotted  in  Fig.  7b.  As  the  amount  of  defocus  decreases  the  standard  deviation  increases 
significantly  due  to  the  narrowing  of  the  PSF.  As  the  blur  spot  narrows  less  of  the  shape  of  the  PSF  can  be 
discerned  from  the  star  images  affecting  the  accuracy  of  the  estimates  of  Z4  .  The  sample  standard  deviation,  os ,  of 
the  joint  parameter  estimate  also  plotted  with  their  associated  CRLB  <y  in  Fig.  7c.  The  CRLB  cr  is  not  achieved  by 
<js  but,  the  standard  deviation  is  below  a  wave  until  the  blur  spot  becomes  too  small. 
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Fig.  7  (a)  The  estimated  defocus  parameter  is  determined  from  simulated  star  data  with  no  noise  present.  The  blue 
X  marks  the  single  parameter  estimate  for  defocus  without  an  atmosphere  present  in  the  simulated  star  data.  The 
green  circles  are  the  single  parameter  estimates  of  Z4  where  r0  =  8 cm  in  the  simulated  star  data.  The  red  boxes  are 


the  joint  parameter  estimate  of  Z4  where  r0  =  8 cm  in  the  simulated  star  data,  (b)  The  joint  parameter  estimates  from 

the  simulated  stars  with  shot  noise  E^Z4  J  and  <js  are  represented  as  blue  dots  and  plotted  as  a  function  of  the 

defocus.  The  FWHM  of  the  PSFs  are  plotted  with  the  green  asterisks  as  a  function  of  defocus,  (c)  The  joint 
parameter  estimates  crs  are  plotted  as  blue  dots  and  the  CRLB  a  as  red  stars  as  a  function  of  star  defocus. 
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6.  CONCLUSIONS 


The  fact  that  the  atmosphere  can  affect  estimation  of  the  coefficient  for  the  defocus  polynomial  demonstrates  the 
need  to  account  for  the  effects  atmosphere  seeing  in  order  to  accurately  estimate  SSTs  aberrations  from  star  images. 
Joint  estimation  of  the  atmospheric  seeing  parameter  for  the  long  exposure  atmosphere  (» 10ms)  and  the  coefficient 

for  the  defocus  polynomial  demonstrates  an  effective  way  to  account  for  the  atmosphere  when  estimating 
aberrations  from  imagery  data.  Results  of  this  work  are  promising  in  that  there  is  no  indication  thus  far  that  joint 
estimation  of  the  Zernike  coefficients  for  the  static  aberrations  in  SST  can’t  be  determined  and/or  alleviated. 

7.  FUTURE  WORK 

The  next  steps  in  this  research  will  be  to  expand  the  joint  estimator  to  include  the  other  aberrations  present  in  SST 
data.  In  addition,  laboratory  experiments  can  be  conducted  by  imaging  point  sources  in  a  diffraction  limited 
telescope  then  defocusing  the  telescope  and  adding  atmospheric  turbulence  to  further  show  that  the  algorithms 
developed  for  joint  estimations  of  Z4  work.  Ultimate  validation  of  the  work  will  be  to  sharpen  the  SST  camera 
images  by  reducing  the  aberrations  in  the  telescope  using  the  basic  techniques  developed  herein. 
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